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ABSTRACT 

We have made high-resolution, high-sensitivity dynamic spectra of a sample of strong pulsars at 
430 MHz with the Arecibo radiotelescope. For 4 pulsars we find faint but sharply delineated features 
in the secondary spectra. These are examples of the previously observed "crisscross" or "multiple drift 
slope" phenomenon presumed to be due to multiple imaging of the pulsar by the interstellar medium. 
The unprecedented resolution and dynamic range of our observations allow a deeper level of analysis. 
Distances to the dominant scattering screen along the line of sight are determined and are shown to 
agree well with those inferred from other scintillation phenomena. Multiple imaging of the pulsar by 
the ISM is required. A compact central image surrounded by a faint scattering halo, roughly circularly 
symmetric, is consistent with the data. Scattering from filaments may also be consistent. The angular 
extent of the scattering material parallel to the direction of the pulsar velocity is ~ 5 mas, corresponding 
to a linear extent of ~ 2 AU. Further observations of these features should allow better discrimination 
between models and an identification of the scattering structures. 

Subject headings: interstellar: matter — pulsars — radio sources: general 



1. INTRODUCTION 



Scattering by density inhomogeneities in the ionized in- 
terstellar medium causes multiple rays to interfere, pro- 
ducing frequency structure in th e spectra of spatially 
coherent sources such as pulsars Scheuer (1968), Rick-| 



ett (1990). The relatively high transverse velocities of 
most pulsars cause the spectra to change significantly on 
timescales of minutes to hours. When the spectrum is 
monitored over time, the resulting two-dimensional array 
or "dynamic spectrum" is often dominated by a random 
pattern of scintillation maxima with characteristic band- 



width and timescale Cordes, Weisberg & Boriakoff (1985) 



On occasion, however, more organized features are present. 
The best studied of these are periodic fringe patterns in 
the dynamic spectra due to the interference between two or 



ing from multiple images on the sky 


Cordes & Wolszczar 


(1986), 


Wolszczan & Cordes (1987)|, 


Gupta, Bhat, & Rac 


(1999), 


Rickett, Lyne & Gupta (1997), hereafter RLG. 



A less well-studied pattern in pulsar dynamic spectra 
consists of a crisscross of intensity maxima that are vari- 
ably spaced perpendicular to the ridge lines of the max- 
ima. Cordes and Wolszczan (1986) drew attention to these 
"multiple drift rates" in dynamic spectra, building on ear- 
lier work (Hewish 1980; Roberts & Abies 1982; Hewish, 
Wolszczan, & Graham 1985; Cordes, Weisberg, & Bori- 
akoff 1985). Numerous examples of the phenomenon are 
evident in the dynamic spectra presented by Gupta, Rick- 
ett, & Lyne (1994) in their multi-year study of 8 pulsars, 
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but the autocorrelation analysis that they employed is not 
well suited to its full elucidation. The remarkable episode 
of "fringing" in the dynamic spectra of PSR B0834+06 
that RLG explore is also accompanied at all epochs by 
this phenomenon. 

The central result of our study is to show that this criss- 
cross pattern has a much clearer representation in trans- 
form space, where it is seen to be a high-Q phenomenon 
with a characteristic parabolic curvature. We discuss the 
constraints that the presence of this pattern places on the 
image geometry and speculate on the physical structures 
that might give rise to this scattering geometry. 

2. OBSERVATIONS AND DATA ANALYSIS 

Our observations were made at the Arecibo Observa- 
tory 3 at 430 MHz. Spectra across a 10 MHz band centered 
on this frequency were produced by a Fourier-transform- 
based spectrometer (AOFTM 4 ). RCP and LCP signals 
were summed in hardware, and spectra with 1024 fre- 
quency channels (i.e. 10 kHz frequency resolution) were 
produced every rs 1.6 ms. During later analysis the spec- 
tra were synchronously averaged at the pulsar period , cal- 

and 
axes 
We used 



culated using TEMPO [Taylor fc Weisberg (1989) 
summed over 10 s time intervals, creating a cube with 
of time, frequency, and pulsar rotational phase. 
512 phase bins across the pulsar period. Dynamic spectra 
were created by producing an ON-pulse spectrum from a 
window of phase bins centered on the pulse and subtract- 
ing an OFF-pulse spectrum to remove the baseline. 
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We observed 19 pulsars during 1999 January and 11 pul- 
sars during 2000 January. A typical observation lasted one 
hour; we made 46 and 43 separate observations at the two 
observing epochs, respectively. Following several previous 
studies (e.g. Cordes & Wolszczan 1986; RLG) we produced 
two-dimensional power spectra from the dynamic spectra. 
These "secondary spectra" form the basis for most of the 
further analysis reported here. We normalized each sec- 
ondary spectrum by its maximum and converted the rela- 
tive power levels into a decibel scale. 

Figure [l] shows dynamic and secondary spectra for 
PSR B0834+06. The crisscross pattern in the dynamic 
spectrum shows up in the secondary spectrum as two 
parabolic features, curving away from the conjugate time 
axis, that we call "arcs". These arcs are remarkable for 
their narrowness and extent. In Figure |^ we present 
secondary spectra, averaged over multiple days to im- 
prove dynamic range. PSR B0823+26 and PSR B0834+06 
showed a persistent arc pattern over both the 1999 and 
2000 observing runs. In 2000, PSR B0919+06 showed arcs 
that were relatively broad and faint and, hence, showed up 
only upon averaging of ~ 6 secondary spectra. However, 
we did not see arcs in two high-quality 1999 observations. 
PSR B1133+16 showed arcs during most of the 1999 and 
2000 observations, but several high-quality 2000 observa- 
tions failed to show them, even though they were made 
on days adjacent to those on which they were present. 
PSR B1933+16 did not show the arc pattern in either of 
two high-quality observations in 1999 or 2000. Time vari- 
ability of the arcs will be the subject of a separate study. 

3. RESULTS 

We detected arcs in the secondary spectra of pulsars 
B0823+26, B0834+06, B0919+06, and B1133+16. By us- 
ing the relationship between fringe frequencies in the sec- 
ondary spectrum and scattering parameters (e.g. Cordes 
& Wolszczan 1986; RLG; Gupta, Bhat, & Rao 1999), we 
can analyze the arcs more quantitatively. Here we summa- 
rize a more detailed study that will be reported elsewhere 
(Cordes et at 2001). To simplify calculations, these treat- 
ments assume that the scattering takes place in a thin 
screen a distance D s from the pulsar, with the pulsar- 
observer distance D. They show that the interference be- 
tween two points #1,2 on the image plane (i.e. as seen by 
an observer) will lead to two-dimensional fringing in the 
dynamic spectrum with fringe frequencies 



ft = 
f» = 



D 
D 

2c~a 



(62 — 9\) • V e fY,±, 

(D-D s )(e 2 2 -el) 



(i) 

(2) 



Here ft is the conjugate time axis, f v is the conjugate fre- 
quency axis, and the effective velocity (Cordes & Rickett 
1998) is a weighted sum of the velocities of the source, 
screen and observer, 

Vea,± = (1 - D S /D)V P± + (D s /D)V obsX - V screen± (3) 
We can introduce the dimensionless variables 

2cD s 



P : 



8i 



D(D - D s ) 



(4) 



q = (02 - 0i) • V e ft,± 



W eff ,j 



(5) 



where V^ff,_i_ is a two-dimensional unit vector for the ef- 
fective velocity. Since p is quadratic in image angles and 
q is linear in image angles, p will generally be quadratic in 
q. This is essential to explaining the parabolic curvature 
of the arcs. 

If we make the reasonable assumption that the pulsar 
velocity is much greater than either the velocity of the 
Earth or the screen, equations (f|)-(||) can be simplified 
to 



P- 



2cs 



D{1 - s) 
As 



Vp ra (l - s) 



ft > 



(6) 
(7) 



where s = D s / D, and we approximate V e g^x by the mea- 
sured proper motion speed, V pm = fj,D. A parabolic arc 
(p = q 2 ) in the p-q plane then corresponds to a feature in 
the secondary spectrum of curvature 



fu 



D A 2 



1 - s 



(8) 



In this model, the curvature of the arcs is determined solely 
by the screen placement, the pulsar distance, and the pul- 
sar proper motion. Thus, by measuring the curvature of 
the arcs in the secondary spectrum, we can determine a 
value for the screen placement that is consistent with the 
observations and assumed values of D and V pm , as in Ta- 
ble 0. 

Gupta (1995) used diffractive bandwidth and timescale 
measurements to determine Vi ss , the speed of the scin- 
tillation pattern past the observer (Cordes and Rickett 
1998). Assuming that the scattering is dominated by a 
thin screen, one can determine the screen placement for 
which Vi SS is consistent with V'pm. This is determined from 
Gupta's Ri parameter by s iss = R 2 /(l + R\)- Table [l] 
shows that the values of Si ss agree well with those for s a ro 
calculated with Eq. 8. Although the values of Si SS agree 
well with those for s arc , the agreement may be only propor- 
tionate because there are significant uncertainties in both 
the distances and the ISS parameters used by Gupta; also, 
scattering is likely to receive contributions from the dis- 
tributed ISM as well as from a screen. In particular, a re- 
cent study of B0919+06 (Chatterjee et al. 2001) that com- 
pares scintillation velocity estimates, which are epoch de- 
pendent, and the parallax and proper motion determined 
from VLB A observations, finds that a thin screen must 
be combined with a distributed medium such as that of 
the Taylor & Cordes (1993) model to bring agreement be- 
tween the two velocity estimates. At some epochs, the 
screen evidently causes multiple imaging of the pulsar. 

From equation (1) it can be seen that the vertical (ft) 
extent of arc features in the secondary spectrum is related 
to the angular extent of the scattering image along the ef- 
fective velocity vector. Table 1 shows the typical angular 
extent of the filament, x ,typ, and the corresponding linear 
size on the screen, xt yp , inferred from the value of ft for a 
generic model. 
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We analyze the distribution of power in the secondary 
spectrum with a sequence of vertical cuts, as in Figure [| 
The rapid decline of power above the arcs strongly limits 
the image geometries that are consistent with the data. In 
order to exclude power above the arc, points must satisfy 
p > q 2 . Referring to equations this constraint can 

be expressed as 



% - e\ > 



?2x 



#lx) 2 , 



(9) 



where the x components are along the effective velocity 
vector. A single extended image has many image pairs 
that violate this condition and, hence, will produce a large 
amount of power above the arcs, inconsistent with the be- 
havior seen. Hence, a multiple image geometry is required 
to produce the clearly delineated arcs we observe. 

One configuration that satisfies Eq. 9, having no de- 
tectable power above the arcs, consists of a bright point- 
like image at the origin (direct line of sight to the pulsar) 
and a faint secondary image as a filled halo around it. The 
interference between a point at the origin and any other 
image point cannot produce power above the arcs because 
Oi = 0. If the point-like image and halo have flux densities 
proportional to a and b, respectively, then the point-like 
source will produce a point-like feature at the origin of the 
secondary spectrum with power proportional to a 2 (e.g. 
Cordes & Wolszczan 1986). The inter- image interference 
between the point-source and the halo will produce power 
proportional to ab. The self-interference between regions 
of the halo (which will produce power above the arcs) will 
be proportional to b 2 . Since this is smaller by another 
factor of b/a from the inter-image term, then, for small 
b/a, the only detectable power will be a point-source at 
the origin and power below the arcs. 



Although the point-source plus halo model is consistent 
with the observations, the scattering material in the screen 
would need to be inhomogeneous enough that the direct 
image of the pulsar would shine through, but rays slightly 
off-axis would be scattered. In essence, we would need to 
be seeing the pulsar through a break in the interstellar 
"clouds" , consistent with the time-variability of the arcs. 
On the other hand, there is strong evidenc e for scattering 



geometries based on filamentary structures Rickctt (1990) 
However, the self-interference between filaments will lead 
to features which are inconsistent with our data. Any 
viable model based on filamentary features in the image 
plane must somehow hide these self-interference features 
as well as produce the observed parabolic arcs. 

The formalism used here includes only geometric delays, 
whereas differential dispersive delays may be important in 
some cases. The fact that the observed arcs conform to the 
geometric-only formalism implies that at least in some in- 
stances the dispersive effects are small. We have examples 
of secondary spectra that show more complex behavior 
and will analyze them within the more general framework 
in a later paper. An extended medium, as opposed to a 
single screen, also will be explored. 

While there are many possibilities for source geometries, 
it is clear that a multicomponent scattered image is nec- 
essary. In future papers, we will apply this method to 
better understand the structure of the ISM on small spa- 
tial scales. 

We thank T. Joseph Lazio for helpful discussions, the 
operators at Arecibo for help with remote observations, 
and Duncan Lorimer for the development of an excellent 
remote observing interface. We acknowledge support from 
the National Science Foundation: grant AST-9618408 to 
JMC and grant AST-9619493 to DRS. 
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Fig. 1. — A dynamic spectrum (top) and its secondary spectrum for a 30 min observation of PSR B0834+06 on 2000 January 18. The 
flux density is displayed as a linear gray-scale with black indicating highest flux. The Nyquist limit along the vertical axis of the secondary 
spectrum is 3 cycles/min, but the spectrum is noise-like beyond the region shown. The gray-scale is logarithmic (linear in decibels), with the 
white level set 3 dB above the noise and the black level set 5 dB below the maximum array value at the origin, as in RLG. These limits span 
a 48 dB range. The strong vertical line in the center of the plot is due partly to the sidelobe response of the Fourier transform and partly 
to broadband pulse-to-pulse fluctuations that are only moderately quenched by the 10-s averaging interval. A similar horizontal line is due 
primarily to the sidelobe response. 

Fig. 2. — Co-added secondary spectra for 5 pulsars. The logarithmic gray-scale is set as in Figure [l| Only a portion of the conjugate time 
axis is shown in several cases. The range of displayed values is 49, 52, 48, 44, 56, and 39 dB for panels (a)-(f), respectively. Panels (c) and 
(d) illustrate the time variability of the arcs; panel (f) is an example of a strong pulsar that does not display arcs. 

Fig. 3. — Cross-sectional cuts through the secondary .spectrum of Figure 0(b), which was first smoothed by 1 cycle/MHz and 0.06 cycles/min 
(10 X 10 pixels). As shown by the tickmarks in Figure bib), the cuts are taken parallel to the vertical axis at steps of 2 cycles/MHz. The cuts 
are made by reflecting the second quadrant of the secondary spectrum through the origin into the fourth quadrant and then cutting vertically 
along the conjugate time axis. The region near where the cuts cross ft = has been set to a constant value since the sidelobes rise to about 
— 10 dB there. Only the first six cuts are labelled. Salient features are the similar height of the arcs on both sides of the origin, the sharp 
drop-off of power above the arcs (i.e. |/t| > 0.5; note the logarithmic axis), the nearly constant power level below the arcs, and the low value 
of the noise floor. 
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Table 1 
Scattering Screen Location 



PSR 


Distance 




f is, axe 


ft, arc 


5 arc 




@x,typ 


x typ 




(kpc) 


(km/s) 


(MHz- 1 ) 


(mill -1 ) 






(mas) 


(AU) 


B0823+26 


0.38 


196 


26 


0.8 


0.36 


0.40 


5.4 


1.3 


B0834+06 


0.72 


174 


31 


0.6 


0.33 


0.36 


4.0 


1.9 


B0919+06 


1.2 


505 


15 


0.55 


0.59 


0.63 


3.6 


1.8 


B1133+16 


0.27 


475 


11 


1.15 


0.49 


0.46 


5.4 


0.7 



Note. — Values are derived from the Taylor and Cordes (1993) distance model and 
the Harrison et al. (1993) proper motions, except for recent interferometric results for 
PSR B0919+06 (Chatterjee et al. 2001). 
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